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TITLE 

TIMING OFFSET COMPENSATION IN ORTHOGONAL FREQUENCY DIVISION 

MULTIPLEXING SYSTEMS 

BACKGROUND OF THE INVENTION 

5 Field of the Invention 

The invention relates to digital communication systems, 
and more particularly to a robust timing offset compensation 
scheme for Orthogonal Frequency Division Multiplexing (OFDM) 
systems . 

10 Description of the Related Art 

With the rapidly growing demand for cellular, mobile 
radio and other wireless transmission services, there has 
been an increasing interest in exploiting various 
technologies to provide reliable, secure, and efficient 

15 wireless communications. Orthogonal Frequency Division 
Multiplexing (OFDM) is well known as a highly spectral 
efficient transmission scheme capable of dealing with severe 
channel impairment encountered in a mobile environment. 
OFDM was previously adopted for wireless local area network 

20 (WLAN) applications as part of the IEEE 802.11a standard in 
the 5 GHz frequency band. Furthermore, the IEEE 802. llg 
standard approved in June of 2 003 also adopted OFDM as a 
mandatory part for a further high-speed physical layer (PHY) 
extension to the 8 02.11b standard in the 2.4 GHz band. 

25 The basic idea of OFDM is to divide the available 

spectrum into several sub- channels (subcarriers) . By making 
all sub-channels narrowband, they experience almost flat 
fading, which makes equalization very simple. In order to 
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obtain a high spectral efficiency, the frequency responses 
of the sub-channels are overlapping and orthogonal. This 
orthogonality can be completely maintained by introducing a 
guard interval, even though the signal passes through a 
5 time-dispersive channel. A guard interval (GI) is a copy of 
the last part of an OFDM symbol which is pre-appended to the 
transmitted symbol. This plays a decisive role in avoiding 
inter-symbol and inter-carrier interference. 

OFDM can largely eliminate the effects of inter- symbol 

10 interference (ISI) for high-speed transmission in highly 
dispersive channels by separating a single high speed bit 
stream into a multiplicity of much lower speed bit streams 
each modulating a different subcarrier. However, OFDM is 
known to be vulnerable to synchronization errors due to the 

15 narrow spacing between subcarriers . The most important 
difficulty when implementing OFDM systems is that of 
achieving timing, phase and frequency synchronization 
between the transmitter and the receiver. In general, 
mismatch between transmitter and receiver oscillators 

20 contributes a non-zero carrier frequency offset in a 
received OFDM signal. Transient behavior of the frequency 
synthesizer is another source of the frequency offset. OFDM 
signals are very susceptible to the frequency offset which 
causes a loss of orthogonality between the OFDM subcarriers 

25 and results in inter-carrier interference (ICI) and bit 
error rate (BER) deterioration of the receiver. Similarly, 
phase noise arising from oscillators also introduces ICI. 
In addition, both frequency offset and phase noise cause 
phase variation so that phase tracking is required for 

30 coherent detection. Unlike the frequency offset and phase 



2 



Client's ref.: P40US 

Our ref.: 07 5 1 - 1 0 2 8 OUS / f ina 1 /M . F . L in/S teve 

noise, as stated above, timing errors may incur inter- symbol 
interference (ISI) in addition to ICI . If the exact timing 
of the beginning of each symbol is not known, the receiver 
cannot be reliably removed the guard interval and correctly 
5 acquire individual symbols before computing the Fast Fourier 
Transform (FFT) of their samples. In this case, inter- 
symbol interference occurs. Moreover, even a small time- 
domain misalignment of the FFT window results in an evolving 
phase shift in the frequency-domain symbols, leading to BER 

10 degradation. Yet another concern issue is the difference 
between the sampling rate of the receiver and that of the 
transmitter. This sampling rate offset results in a 
rotation of the 2 m - ary constellation from symbol to symbol. 
Many techniques dealing with the frequency offset 

15 estimation have been previously proposed for OFDM systems. 
Nevertheless, fewer works lend themselves readily to timing 
and phase compensation suitable for integrated circuit 
implementation. In order to achieve rapid acquisition and 
accurate tracking, there is a need to particularly address 

20 the timing offset in the lock transient of local 
oscillators . 

SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a 
timing estimation apparatus for mult i -carrier systems such 
25 as OFDM receivers. 

It is another object of the present invention to 
provide a timing tracking apparatus for rapidly acquiring 
and continuously tracking the timing of OFDM signals. 
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It is yet another object of the present invention to 
provide a robust timing offset compensation scheme for OFDM 
systems, which is well-suited to integrated circuit 
implementation . 

5 The present invention is generally directed to a timing 

offset compensation apparatus for multi-carrier systems. 
According to one aspect of the invention, a timing error 
estimation apparatus for multi-carrier systems is disclosed. 
The timing error estimation apparatus includes a timing 

10 offset compensator and a timing error estimator. The timing 
offset compensator receives a current symbol in the 
frequency domain after taking an Af-point Discrete Fourier 
Transform (DFT) and compensates the current symbol for the 
effect of timing offset with a timing offset prediction 

15 value. The timing error estimator is coupled to the timing 
offset compensator to take a timing compensated version of 
the current symbol on pilot subcarrier locations; it can 
calculate a timing error value based on a function of a 
phase tracking value, a channel response of each pilot 

20 subcarrier, transmitted data on each pilot subcarrier, and 
the timing compensated version of the current symbol on said 
pilot subcarrier locations. 

According to another aspect of the invention, a timing 
tracking apparatus for multi-carrier systems is proposed. 

25 The timing tracking apparatus is constituted by a parameter 
table, an «th-order tracking loop and a timing synchronizer. 
The parameter table is configured to store a number of loop 
parameters. The «th-order tracking loop is adapted to 
compute a timing offset prediction value for a next symbol 

30 based on a timing error value of a current symbol, the 
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timing offset prediction value of the current symbol and the 
loop parameters that are retrieved from the parameter table 
for the current symbol. According to the timing offset 
prediction value of the next symbol, the timing synchronizer 
5 generates a shift amount of the DFT window for the next 
symbol, in which the shift amount is equal to zero if the 
timing offset prediction value of the next symbol is within 
a predetermined range. Moreover, the timing synchronizer 
applies the DFT window shift amount to align the DFT window 

10 and further starts an inhibit interval if the shift amount 
is not equal to zero. Upon completion of the inhibit 
interval, the DFT window shift amount is further subtracted 
from the timing offset prediction value. 

In a preferred embodiment of the invention, a timing 

15 offset compensation apparatus is composed of a timing offset 
compensator, a timing error estimator and a timing tracking 
unit. The timing offset compensator receives a current 
symbol in the frequency domain after taking an //-point 
Discrete Fourier Transform (DFT) and compensates the current 

20 symbol for the effect of timing offset with a timing offset 
prediction value. The timing error estimator is adapted to 
take a timing compensated version of the current symbol on 
pilot subcarrier locations; it calculates a timing error 
value for the current symbol based on a function of a phase 

25 tracking value, a channel response of each pilot subcarrier, 
transmitted data on each pilot subcarrier, and the timing 
compensated version of the current symbol on the pilot 
subcarrier locations. The timing tracking unit receives the 
timing error value of the current symbol to generate the 
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timing offset prediction value and a shift amount of the DFT 
window for a next symbol . 

DESCRIPTION OF THE DRAWINGS 

The present invention will be described by way of 
5 exemplary embodiments, but not limitations, illustrated in 
the accompanying drawings in which like references denote 
similar elements, and in which: 

FIG. 1 is a flowchart illustrating primary steps of 
dealing with feedback delay from the timing offset tracking 
10 to the FFT (DFT) window alignment; 

FIG. 2 is a block diagram of a preferred embodiment 
according to the invention; and 

FIG. 3 is a detailed block diagram of a timing tracking 
unit according to the invention. 

15 DETAILED DESCRIPTION OF THE INVENTION 

The present invention will now be described in the 
context of the use of OFDM for communication, although the 
present invention is not limited to OFDM. The present 
invention is also described with reference to a wireless 
20 communication system that conforms to the IEEE 802.11a/g 
standard. According to the invention, the communication 
system need not be wireless and the conformant 802.11a/g 
transceiver referred to herein is merely an exemplary 
illustration . 

25 In a conformant 802.11a/g system, an OFDM symbol is 

modulated onto a number of subcarriers by applying an Ap- 
point Inverse Fast Fourier Transform (FFT) with N=64. At 
the receiving end, the demodulation can be accomplished by 
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means of an N- point FFT. Ignoring the ICI noise term, the 
output of FFT is described by 

R,, k =H k X i>k e»>e~ J * * +N iJc (1) 

and 

5 ^ = n-/+^ 0/ n = 27f £ T f 

where 

i denotes the OFDM symbol index, 
k denotes the subcarrier index, 

X ik is complex-valued data transmitted on subcarrier k 
10 of the ith OFDM symbol, 

H k is the channel response of subcarrier k, 

N iJc is the frequency- domain representation of additive 

white Gaussian noise (AWGN) at the Arth subcarrier 

during the ith OFDM symbol , 
15 T is the FFT period, 

At i is the timing offset with respect to the FFT period 

during the ith OFDM symbol, 

V is the symbol interval (including the GI duration) , 
f e is the residual frequency error due to the initial 
20 frequency offset estimation, 

0 O is the initial phase drift, and 

0,. is the equivalent phase offset of the ith OFDM 
symbol . 

Note that the noise term, N ik , is iid (identically, 
25 independently distributed) over different symbols and 
subcarriers. From the above signal model, the phase and 
timing offsets both result in additional phase shifts. It 
is easily seen that the presence of timing offset causes a 
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phase rotation proportional to the subcarrier index k and 
the timing offset At i . 

Denoting the pilot subcarrier index by p m , the 
conditional joint probability density function of those 
5 pilot tones of the ith symbol given fa , At, , H k and X ik , for 
k~ p m , is given by: 



r 



1 



exp 



1 % 



R -H X e J(l>i e T 



2\ 



where 



10 



2 

<j n denotes the noise variance per subcarrier and 



denotes the ith received OFDM symbol on pilot 



subcarrier location p m , and 

N SP is the number of pilot subcarriers. 
Therefore, a pilot-aided estimation is derived from this 
15 function using the maximum- likelihood (ML) principle as 
follows : 

A?, = arg max gfa^ , R ip% , ■ . . , R.^ | fa , At, , H k , Jf a ) 



= arg max 
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L(t)=Re 



m=\ 



V J 



Note that superscript herein denotes complex conjugation. 

To find the maximum of L{t) , the equation of = !,'(/)= 0 

dt 

5 must be solved first. However, the derivative L'(t) can be 
shown to have the form: 



L'(t)=Re 



n j2n— » 



= ^Re 
T 



f 



j Z p m R iaf 



PJ 
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Im 



2 Pm R U k 



-/2/r 
n p m *i*P m e e 

V 



H X e Jipi e T 



Although it is hard to find the close form solution to the 
equation L'(t) = O f there are many numerical methods suitable 
10 for finding an approximation to the solution. For example, 
Newton's method is an efficient iterative method for 
estimating solutions to equations. Using Newton's method, 
an approximate solution t n+x to the equation L'(t) = 0 can be 
obtained from a previous estimate t n as follows: 



where the second derivative of L(t) is given by: 



(2) 
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L"(t) = -f Im 



No, 
m=l 



J y Pm ' pn > 



7" 



Im 



Re 



m=\ 



£ p m R,, f 

m=l 



-72^ 



Note that the second derivative L n (t) is in the denominator 
of equation (3) . This may lead to strong noise enhancement 
at low SNR. To mitigate the noise effect, the major noise 
5 source in L"(t) namely R ( should be addressed. From 
equation (1) , /? f . can be approximated by dropping the noise 



term N ik 



-JL71— L 



Substitution this approximation into L"(t) yields: 



10 



Due to \ X i, Pm \ = 1 # the second derivative L n (t) can be further 
reduced to: 

Here equation (2) is normalized by a sampling period of 
15 T s = T/N : 



/ n+1 _ t„ 1 L'{t„) 
T s T s T s L"(t„) 

This can be rewritten as: 
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i r(r,,rj 

where 

r Ei is a measure of timing offset of the ith symbol, and 
r Pi is a timing offset prediction value of the ith 
5 symbol . 

According to the invention, a second-order tracking loop is 
employed to follow the transmitter's behavior during the 
lock transient and continuously track the timing of OFDM 
signals. The tracking loop is modeled with a set of 
10 recursive equations by: 

T T,i = T P,i + MrJ T £,i , % 

(4) 

and 
where 

15 r Ti and v Ti denote timing and period offset tracking 

values of the ith symbol, respectively, 

ju Ti and fi vi denote loop parameters of the ith symbol 
for t tj and v Ti , respectively, 

T Fi+x is the timing offset prediction value of symbol 

2 0 1+1, 

v T4 _ x is the period offset tracking value of symbol i-l , 
and r ei , a timing error value of the ith symbol, is given by: 

T £j ~ T E,i ~ T P,i 

Note that variables or parameters with subscript i refer to 
25 the current symbol and those subscripted with i+l refer to 
the next symbol. On the other hand, variables or parameters 
with subscript i-l indicate that they are associated with 
the preceding symbol. The index of i=0 represents the first 
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symbol to be tracked in the tracking loop, which is the 
SIGNAL symbol in IEEE 802.11a. Initial values of recursive 
equations (4) and (5) are indexed with i = -l. According to 
the invention, the initial values are given by: 

t t-\ = 0 

J c 



and 



r p,o ~" t t-\ + v t-\ ~ V T-\ 



where f d is an estimate of frequency offset, f c is the 
nominal carrier frequency, and N' = T'/T S . Furthermore, the 
10 tracking of period offset should skip the first symbol. 
Thus , fd v 0 = 0 . 

The phase shift caused by timing offset A*,, can be 
compensated for by means of multiplying the FFT output of 

the /th symbol, R ik , by a compensation factor of e N . Let 

15 r i -bd i jT s ; timing offset At i in equation (1) is now 
normalized to the sample in respect of the receiver's FFT 
window. If x i exceeds one sample, FFT window alignment will 
be required. According to invention, the shift amount for 
the FFT window alignment is given by: 



2 o *dj + l = 



+ 1* Sfljy+i >r v 
0, otherwise 



where T K and F A are thresholds for decision. The shift 
amount of A rf l+I = +1 means that the FFT window should be 
delayed by one sample in time, while A rf)I . +1 =-l means that the 
FFT window should be advanced by one sample. To address the 
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unsteady timing shifts, a hysteresis range is introduced to 
the above criterion. Therefore, 

r K -r x =i + r H 

where Y H namely the hysteresis range is a positive real 
number. In addition, it is necessary to deal with feedback 
delay from the timing offset tracking to the FFT window 
alignment. FIG. 1 illustrates a way to handle this issue. 
Briefly, A rf>1+1 is applied to align the FFT window for the 
next symbol and an inhibit interval will be started if A d M 
is not equal to zero. Upon completion of the inhibit 



interval , 



V,/+i 



is provided for further subtraction from 



'/> i + l 



In order to prevent the FFT window from additional 



adjustments, however, A di+l is set to 0 during the inhibit 
interval. It should be contemplated now how to estimate 
15 timing error value r £i for the ith symbol. On the basis of 
equation (3), r ei can be calculated from 



1 L'(t pj T s ) 
EJ T s L»{t pj T s ) 



In. 

Im 

T 



2 P m R i, f 
m=l 



H X. e j(j>i e N 



V 



J 



N 



lm< e 



m=\ J 



(2n 



(6) 



In 



n sp 2 I I 2 
Pm\ H p m \ 

m=\ 



where R[ Pm , a timing compensated version of the ith received 
OFDM symbol on pilot subcarrier location p m , is given by: 

Pm T PJ 



R U=* ( ,p m e " 



(7) 
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While the FFT is mentioned in the above discussion, it 
should be clear to those skilled in the art that the 
Discrete Fourier Transform (DFT) is also applicable to the 
present invention since the FFT is an efficient scheme for 
5 computing the DFT. Therefore, DFT and FFT are herein 
interchangeable terms according to the principles of the 
invention . 

In light of the foregoing description, the present 
invention will now be explained from an embodiment of FIG. 

10 2 . As depicted, an FFT (DFT) processor 250 accepts a 
baseband signal at its input and applies an N- point FFT 
(DFT) to transform the baseband signal from the time domain 
to the frequency domain. The FFT (DFT) output R ik is then 
fed to a timing offset compensation apparatus of the 

15 invention. The timing offset compensation apparatus 200 
includes a timing offset compensator 210, a timing error 
estimator 220 and a timing tracking unit 230. The timing 
offset compensator 210 receives a current symbol, R ik , from 
the FFT (DFT) processor 2 50 and compensates R ih for the 

20 effect of timing offset with a timing offset prediction 
value, r Pi . After that, the timing error estimator 22 0 
takes R'. , a timing compensated version of the current 
symbol on pilot subcarrier locations p m where m=l,..., N SP . 
In addition, the timing error estimator 22 0 also takes known 

25 data X iPm that are transmitted on pilot subcarriers of the 
current symbol. A channel estimator 260 is used to estimate 
the channel response of each subcarrier, H k . The channel 
response of pilot subcarriers, H Pm and /w=l,..., N sp , is then 
provided to the timing error estimator 220. Moreover, a 

30 phase tracking loop 2 70 is configured to generate a phase 
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tracking value , <f> Ti , for the current symbol . When 
generated, <j> Ti is fed to the timing error estimator 220. As 
a result, the timing error estimator 22 0 calculates a timing 
error value r eJ for the current symbol with equation (6) in 
5 which fa is replaced by </> Ti so equation (6) becomes a 
function of </> TJ , H Pm , X i%Pm , and R\ Pm . 

The timing tracking unit 23 0 is designed to generate 
the timing offset prediction value r PM and a shift amount 
of the FFT (DFT) window, A dJ+l , for a next symbol based on a 

10 frequency offset estimate & d from outside and the timing 
error value r eJ . Referring to FIG. 3, a block diagram of 
the timing tracking unit 230 is illustrated. The timing 
tracking unit 230 is constituted by a parameter table 310, 
an rcth-order tracking loop 32 0 and a timing synchronizer 

is 330. The parameter table 310 is configured to store a 
number of loop parameters. The nth-order tracking loop 320 
generates the timing offset prediction value r Pi+l for the 
next symbol based on the timing error value r €i , a prior 
timing offset prediction value r Pi and the loop parameters 

20 ju Ti and p vi that are retrieved from the parameter table 310 
for the current symbol. The nth-order tracking loop 320 
internally computes a timing offset tracking value r Ti and a 
period offset tracking value v Ti for the current symbol 
based on r Pi , r £i , /u Ti and ju vi , and yields r PM by summing 

25 r Ti and v Ti . In one embodiment, the wth-order tracking loop 
32 0 is a second-order tracking loop modeled with the above- 
described equations (4) and (5) . 

Turning now to FIG. 1 where substantial execution steps 
of the timing synchronizer 330 are illustrated. In step 

30 S101, DELAY_CNT is checked to see whether it is equal to 
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zero, in which DELAY_CNT is used to do a countdown of an 
inhibit interval. If DELAY_CNT is zero (which means the 
inhibit interval is not activated) , then the process 
proceeds to step S103 where the timing synchronizer 330 
5 generates the FFT (DFT) window shift amount for the 

next symbol according to t pj+x . In the meantime, the timing 
synchronizer 33 0 keeps A di+l in 8 shift . Note that A rf/+1 is set 
to zero if r Pi+l derived from the tracking loop 32 0 is within 
a predetermined range between T K and F A . In step SI 05, 

10 & shift * namely A d M , is checked to determine whether it is not 
equal to zero. If so, the process proceeds to step 107 
where the timing synchronizer 330 starts the inhibit 
interval and sets DELAY_CNT to a predetermined value of 
TIMING_FB_DLY. Otherwise, the process proceeds to step 111 

is where DELAY_CNT is checked to see whether it is equal to 
one. In step S115, T Pi+x is directly provided as output if 
DELAY_CNT is not equal to one. On the other hand, 8 shift is 
subtracted from r PM in step S113 if DELAY_CNT is equal to 
one. Subsequently, S shift is reset to zero. In step S117, 

20 DELAY_CNT is checked to see whether it is greater than zero. 
If so, DEIiAY_CNT is decreased by one; otherwise, the process 
proceeds back to step S101 for a new symbol. When DELAY_CNT 
is not zero, A d>f+I is set to zero and the process proceeds 
to step Sill. As such, the FFT (DFT) window is prevented 

25 from further adjustments during the inhibit interval. In 
summary, the timing synchronizer 33 0 applies A dJ+{ to align 
the FFT (DFT) window for the next symbol and further starts 
the inhibit interval if A dM is not equal to zero. 
Additionally, S $hift (the effective FFT/DFT window shift 



16 



Client's ref . : P40US 

Our ref.: 0751-10280 US /final /M . F . Lin /Steve 

amount) is subtracted from r Pi+} immediately upon completion 
of the inhibit interval . 

Referring again to FIG. 3, the timing tracking unit 230 
further comprises a D flip-flop 322 that receives r PM from 
5 the timing synchronizer 33 0 and provides as output the 
timing offset prediction value, x Pi , for the current symbol. 
As a result, the timing offset compensator 210 can 
compensate the current symbol using r PJ and provides as 
output a timing compensated version of the current symbol, 

10 R' ik . In this way, the apparatus of the invention really 
addresses the timing offset in the lock transient of local 
oscillators. Accordingly, the present invention provides a 
robust scheme to rapidly acquire and continuously track the 
timing of OFDM signals. 

15 While the invention has been described by way of 

example and in terms of the preferred embodiments, it is to 
be understood that the invention is not limited to the 
disclosed embodiments. To the contrary, it is intended to 
cover various modifications and similar arrangements (as 

20 would be apparent to those skilled in the art) . Therefore, 
the scope of the appended claims should be accorded the 
broadest interpretation so as to encompass all such 
modifications and similar arrangements. 
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